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I.  Introduction 


Since  the  construction  of  the  Large  Blast  Simulator  (LBS)  in  Gramat,  France  [1],  the 
United  States  has  been  actively  engaged  in  developing  a  design  for  a  Large  Blast/Thermal 
Simulator  (LB/TS)  [2-10].  The  primary  purpose  of  this  part  of  the  design  study  was  to 
estimate  the  blast  loading  on  a  candidate  passive  Rarefaction  Wave  Eliminator  (RWE)  design 
for  the  proposed  US  LB/TS.  This  passive  RWE  design  concept  arid  others  were  conceived 
and  tested  at  small  scale  at  the  BRL  [11],  and  scaled  up  to  the  LB/TS  [12].  The  RWE  is 
located  at  the  discharge  end  of  the  LB/TS.  Its  purpose  is  to  prevent,  or  at  least  minimize, 
the  generation  of  upstream-traveling  rarefaction  waves  caused  by  the  interaction  of  the  blast 
wave  with  the  ambient  air  outside  of  the  LB/TS  by  a  combination  of  partially  reflecting  the 
incident  shock  and  by  acting  as  a  valve  to  control  the  flow.  Otherwise,  these  waves  could 
travel  upstream  against  a  subsonic  flow  and  distort  the  remainder  of  the  blast  wave  at  the 
test  station.  A  simplified  schematic  of  the  proposed  US  LB/TS  is  shown  (not  to  scale)  in 
Figure  la,  and  a  cross-section  of  the  LB/TS  is  shown  in  Figure  lb.  This  design  includes 
an  Interior  (upstream)  RWE  (hereinafter  referred  to  as  the  I- RWE)  that  has  a  set  amount 
of  blockage  caused  by  fixed  vertical  bars  equal  to  20  percent  of  the  cross-sectional  area,  of 
the  LB/TS  test  section,  and  a  Primary  RWE  (hereinafter  referred  to  as  the  P-RWE)  with 
adjustable  blockage  settings  from  31.5  percent  to  80  percent  of  the  LB/TS  cross-sectional 
area  at  the  exit  of  the  expansion  tunnel  (see  Figure  2).  A  drawing  of  the  P-RWE  is  shown 
in  Figure  3,  with  cross-sections  of  two  generic  candidate  shapes  for  a  single  vane  shown 
in  Figure  4'.  A  drawing  of  the  I-RWE  is  shown  in  Figure  5.  Also  included  in  this  design 
are  two  sets  of  side  vents  which  can  be  either  completely  closed,  or  opened  to  a  combined 
maximum  of  61.5  percent  of  the  LB/TS  cross-sectional  area  (see  Figure  2).  The  set  of  side 
vents  located  between  the  I-RWE  and  the  P-RWE  is  termed  the  Downstream  Side  Vent 
(hereinafter  referred  to  as  the  DSV),  and  the  set  of  side  vents  located  upstream  from  the  I- 
RWE  is  termed  the  Upstream  Side  Vent  (USV).  The  DSV  is  necessary  to  provide  additional 
venting  area  to  compensate  for  the  unavoidable  blockage  (a  minimum  of  31.5  percent)  of 
the  P-RWE  in  its  fully  opened  position.  The  USV  provides  similar  additional  venting  to 
compensate  for  the  fixed  20  percent  blockage  of  the  I-RWE. 

A  second  purpose  of  this  study  was  to  estimate  the  overpressure  buildup  at  the  end  of 
the  LB/TS  near  the  RWE’s,  particularly  in  a  situation  where  the  combined  open  areas  of  the 
P-RWE,  USV,  and  DSV  are  too  small  for  the  planned  blast  wave.  The  strongest  blast  wave 
currently  planned  for  simulation  in  the  LB/TS  is  a  241.3  kPa  (35  psi)  peak  overpressure 
blast  wave  from  a  600  KT  nuclear  burst  above  ground,  so  this  wave  was  chosen  for  this 
study.  Two  hydrocode  computations  were  run,  one  assuming  “correct”  P-RWE  and  side 
vent  flow  area  settings  equal  to  a  total  of  130  percent  of  the  LB/TS  cross-sectional  area*  and 
ariother.  assuming  an  accidental  “worst  case”  P-RWE  and  side  vent  total  flow  area  setting 
equal  to  20  percent  of  the  LB/TS  crOss-sectional  area.  The  computation  using  the  correct 
P-RWE  settings  was  assigned  probleiri  number  8606.02.  For  simplicity,  this  computation 
will  hereinafter  be  referred  to  as  the  “normal  vent”  computation.  The  computation  using 
the  presumed  worst  case  vent  area  was  assigned  problem  number  8605.27.  For  simplicity, 
this  computation  will  hereinafter  be  referred  to  as  the  “low  vent”  computation.  These 
computations,  are  described  in  the  next  section. 
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II.  The  Hydrocode  Computations 

The  BRL  version  [13-14]  of  the  HULL  [15]  hydrodynamic  computer  code  was  used  for  the 
numerical  propagation  of  the  blast  wave  down  the  LB/TS  and  its  interaction  with  the  RWE’s 
and  side  vent  areas.  This  version  of  the  HULL  hydrocode  solves  the  inviscid  Euler  equations 
by  using  a  modified  Lax-Wendroff,  explicit  time  step,  finite  difference  method.  HULL  is 
actually  a  large  collection  of  cooperative  programs  and  shared  subroutines.  Its  principal 
programs  are  a  grid  generator  KEEL,  the  hydrocode  HULL,  and  a  plotting  utility  PULL. 
This  version  of  HULL  can  be  set  up  to  run  using  either  three-dimensional  (3-D)  or  two- 
dimensional  (2-D)  Cartesian  coordinates,  or  2-D  cylindrical  coordinates  having  rotational 
symmetry.  All  flow  field  cells  in  3-D  Cartesian  coordinates  must  be  either  entirely  filled 
with  hydrodynamic  material  or  designated  as  non-hydrodynamic  and  rigid  (“island”)  cells. 
Flow  field  cells  in  the  2-D  versions  may  also  be  of  either  type,  and  may  additionally  be  half 
hydrodynamic  and  half  island  (“shore”)  [14].  All  cells  in  all  versions  must  have  a  rectangular 
cross-section; 

In  order  to  make  the  hydrocode  computations  more  tractable,  a  geometric  simplification 
was  made.  Although  the  problem  in  the  actual  LB/TS  is  inherently  3-D  (see  Figure  1),  most 
of  the  principal  features  of  the  flow  for  this  problem  of  determining  RWE  blast  loading  can  be 
approximated  well  in  2-D  HULL  with  cylindrical  coordinates.  The  problem  was  transformed 
into  an  axially  symmetric  problem,  with  a  cylindrical  shock  tube  having  a  radius  equal  to 
the  hydraulic  radius  of  the  LB/TS.  The  hydraulic  radius  is  equal  to  the  cross-sectional  area 
of  the  LB/TS  divided  by  its  perimeter.  Flow  and  blockage  areas  for  the  I- RWE  and  P-RWE 
were  transformed  into  spaced  concentric  toroidal  rings  of  either  hydrodynamic  or  island  cells, 
and  the  two  bands  of  discrete  side  vents  which  comprise  the  USV  and  DSV  were  transformed 
into  radial  vent  bands  around  the  entire  circumference  of  the  transformed  LB/TS.  The  side 
wall  of  the  LB/TS,  except  for  the  open  bands  representing  the  side  vents,  was  simulated  by  a 
rigid,  impermeable  shell  made  up  of  island  cells.  The  drivers  were  not  simulated,  but  rather 
the  LAMB  [16]  code,  a  set  of  empirical  fits  based  in  part  on  the  1  KT  Standard  [17],  was  used 
as  the  upstream  input  boundary  condition  to  feed  the  241.3  kPa,  600  KT,  decaying  blast 
wave  into  the  HULL  hydrocode  grid  as  time  progressed.  A  similar  feed-in  of  a  LAMB  blast 
wave  through  the  same  type  of  transmissive  boundary  is  discussed  in  a  previous  study  [18] 
on  shock  tube  blockage  effects.  The  loading  on  the  RWE’s  is  produced  not  only  by  the 
direct  effect  of  the  interaction  of  the  blast  wave,  but  is  also  enhanced  by  the  blockage  effects 
documented  in  previous  work  [18-22]. 

The  low  vent  hydrocode  computation  was  for  the  worst  case  situation  where  both  the 
USV  and  DSV  were  closed  and  the  P-RWE  was  only  open  to  20  percent  of  the.  LB/TS 
cross-sectional  area  (i.e.,  80  percent  closed).  The  normal  vent  hydrocode  computation  was 
for  the  situation  where  the  P-RWE  and  the  side  vents  were  set  at  their  design  positions 
for  this  blast  wave:  68.5  percent  open  for  the  P-RWE,  16.6  percent  of  the  test  section 
area  for  the  DSV,  and  44.9  percent  of  the  test  section  area  for  the  USV.  These  openings 
all  represent  maximum  possible  openings  under  this  design  concept.  The  I-RWE  flow  area 
was  fixed  at  80  percent  of  the  LB/TS  cross-sectional  area  (i.e.,  20  percent  closed)  for  both 
computations.  Both  computations  used  the  same  basic  grid  design,  including  having  identical 
flow  field  computational  cell  sizes,  with  the  only  variations  being  those  necessary  to  simulate 
the  differences  in  side  vent  areas  and  P-RWE  settings.  These  differences  were  a  matter  of 
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defining  some  cells  in  the  LB/TS  side  wall  as  either  hydrodynamic  cells  for  the  normal  vent 
case  through  which  material  could  flow,  or  island  cells  for  the  low  vent  case  which  created  a 
barrier  to  the  flow.  A  large  free  field  was  defined  outside  of  the  LB/TS  so  that  outflow  from 
the  LB/TS  would  not  be  artificially  interfered  with  by  the  computational  boundaries  of  the 
grid. 

It  should  be  noted  that  the  exact  dimensions  and  locations  used  in  these  computations 
are  not  necessarily  design  values  for  the  LB/TS,  but  in  some  cases  are  the  results  of  the 
geometric  simplifications  introduced  into  the  problem  to  make  a  more  tractable  computation 
using  a  cylindrically  symmetric  geometry  in  the  HULL  hydrocode.  The  essentials  of  the  grids 
in  the  axial  direction,  a  total  of  458  cells,  are  listed  below: 

1.  A  210  cell  region  defined  as  Set  A,  from  the  start  of  the  grid  (axially)  at  0.0  m  to 
51.7120  m,  which  includes  the  test  station  at  20.0000  m.  (Note  that  this  zero  point  is 
not  the  true  start  of  the  LB/TS  expansion  section,  but  is  used  simply  as  a  matter  of 
convenience  for  these  computations.  The  quotation  of  four  significant  figures  beyond 
the  decimal  point  is  done  to  give  an  indication  of  the  sizing  of  the  flow  field  cells.)  The 
normal  wall  thickness  in  the  LB/TS  ends  at  5L7120  m,  which  is  the  beginning  of  the 
thick  wall  section  at  the  RWE  end  of  the  expansion  section.  The  current  design  for  the 
LB/TS  calls  for  the  last  several  meters  of  the  expansion  section  to  be  approximately 
1.22  m  (4.0  feet)  thick.  This  section  serves  in  an  additional  capacity  as  a  carrier  of 
anchoring  cables  for  the  P-RWE.  The  thickness  of  the  LB/TS  wall  is  not  directly  relevant 
to  these  computations,.  It  is  modeled  as  0.862  m  thick  throughout  its  entire  length. 

2.  A  69  cell  region,  Set  B,  from  51.7120  m  to  62.4659  m,  which  includes  the  thick  wall 
section  upstream  from  the  USV.  The  HULL  hydrocode  computation  is  started  with  the 
shock  front  located  in  this  Set  at  61.9575  m  for  both  computations. 

3.  A  16  cell  region,  Set  C,  from  62.4659  m  to  64.0929  m,  whose  intersect  with  Set  L  (see 
the  radial  grid  description  below)  simulates  the  open  span  in  the  LB/TS  side  wall  for 
the  USV  in  the  normal  vent  problem.  The  open  area  is  equal  to  44.9  percent  of  the 
LB/TS  cross-sectional  area.  The  intersect  of  this  Set  with  Set  L  is  a  solid  wall  for  the 
low  vent  problem. 

4.  A  nine  cell  region,  Set  D,  from  64.0929  m  to  65.0000  m,  which  simulates  a  solid  side 
wall  region  between  the  USV  and  the  I-RWE. 

5.  A  six  cell  region,  Set  E,  from  65.0000  m  to  65.6600  m,  which  simulates  the  I-RWE, 
which  blocks  20  percent  of  the  LB/TS  cross-sectional  area  for  both  computations.  This 
blockage  is  simulated  by  a  grouping  of  five  concentric  toroidal  rings  of  island  cells  within 
the  intersect  of  this  Set  and  Set  K. 

6.  A  20  cell  region,  Set  F,  from  65.6600  m  to  67.6993  m,  which  simulates  the  region  between 
the  I-RWE  and  the  DSV. 

7.  A  six  cell  region,  Set  G,  from  67.6993  m  to  68.3008  m,  whose  intersect  with  Set  L 
simulates  the  open  area  in  the  LB/TS  side  wall  for  the  DSV  in  the  normal  vent  problem. 
The  open  area  is  equal  to  16.6  percent  of  the  LB/TS  cross-sectional  area.  The  intersect 
of  this  Set  with  Set  L  is  a  solid  wall  for  the  low  vent  problem. 
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8.  A  17  cell  region,  Set  H,  from  68.3008  m  to  70.0000  m,  which  simulates  the  region 
between  the  DSV  and  the  P-RWE  at  the  end  of  the  LB/TS. 

9.  An  eight  cell  region,  Set  I,  from  70.0000  m  to  70.8440  m,  which  simulates  the  P-RWE. 
This  is  set  to  an  opening  equal  to  20  percent  of  the  LB/TS  cross-sectional  area  for  the 
low  vent  problem  (i.e.,  80  percent  closed).  It  is  set  to  an  opening  equal  to  68.5  percent  of 
the  LB/TS  cross-sectional  .area  for  the  normal  vent  problem  (i.e.,  31.5  percent  closed.) 
These  different  blockages  are  simulated  by  different  choices  of  groupings  of  toroidal  rings 
of  island  cells  within  the  intersect  of  this  Set  and  Set  K. 

10.  A  97  cell  region,  Set  J,  from  70.8440  m  to  102.389  m,  which  simulates  the  distance  axially 
from  the  end  of  the  P-RWE  to  the  far  downstream  boundary  in  the  axial  direction. 

The  essentials  of  the  grids  in  the  radial  direction,  a  total  of  144  cells,  are  listed  below: 

1.  A  60  cell  region,  Set  K,  from  0.0  m  to  7.24715  m,  which  simulates  the  interior  of  the 
.  LB/TS  from  the  origin  to  the  LB/TS  side  wall.  (Note  that  7.24715  m  is  the  hydraulic 

radius  of  the  actual  LB/TS.) 

2.  A  six  cell  region',  Set  L,  from  7.24715  m  to  8.10925  m,  which  simulates  the  LB/TS  side 
’  •  wall,  typically  with  island  cells  except  where  side  vents  are  simulated.  This  thickness 

was  chosen  to  approximate  the  cross-sectional  area  of  the  side  wall  if  it  were  uniformly 
distributed  in  a  complete  circular  shell.  This  was  an  arbitrary  decision.  To  keep  the 
grid  simple,  the  side  wall  was  simulated  to  be  of  the  same  thickness  throughout  its 
length;  any  effect  on  the  loadings  of  interest  here  because  of  this  simplification  would 
be  negligible. 

3. '  A  78  cell  region,  Set  M,  from  8.10925  m  to  40.1679  m,  which  simulates  the  region 

between  the  outer  surface  of  the  LB/TS  side  wall  and  the  outer  radial  boundary  of  the 
computational  flow  field. 


III.  Loading  on  the  Rarefaction  Wave  Eliminators 

1.  Primary  RWE  . 

a,.  Overpressure  on  the  Primary.  RWE 

The  overpressure  on  the  P-RWE  for  the  low  vent  problem  is  shown  in  Figure  6a.  The 
values  for  average  overpressure  on  the  front  and  back  faces  are  computed  using  only  those  flow 
field  cells  immediately  adjacent  to  the  island  cells  simulating  the  P-RWE  itself.  These  values 
are  then  used  to  obtain  the  net  axial  average  pressure  acting  on  the  P-RWE,  subtracting 
the  overpressure  on  the  back  from  that  on  the  front.  The  net.  overpressure  curve  shows  a 
rapid  rise  to  a  peak  average  overpressure, .of  658.9  kPa  (95.6  psi)  at  85.4  ms,  not  too  far 
from  the  peak  theoretical  overpressure  of  850.2  kPa  (123.3  psi)  for  normal  reflection  from 
a  continuous  wall.  (Values  qouted  here  and  in  subsequent  discussions  are  from  tabulated 
results.)  There  is  then  a  subsequent  rise  to  704.4  kPa  (102.2  psi)  at  108.4  ms  because  of  the 
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return  of  a  partial  reflection  from  the  back  face  of  the  I-RWE,  and  a  slow  reduction  in  net 
pressure  from  that  point.  The  time  used  for  this  plot,  and  all  subsequent  plots,  is  referenced 
to  t  =  0.0  at  the  incident  shock  arrival  at  the  test  station,  located  in  both  computations  at 
the  20.0  m  axial  position  as  defined  for  purposes  of  these  computations. 

The  overpressure  on  the  P-RWE  for  the  normal  vent  problem  is  shown  in  Figure  6b. 
The  time  of  arrival  of  the  shock  is  essentially  the  same  as  for  the  low  vent  problem,  but  the 
peak  net  overpressure  on  the  P-RWE  is  476.8  kPa  (69.2  psi),  28  percent  below  the  first  peak 
for  low  venting.  This  lower  peak  is  also  relieved  very  rapidly,  with  no  second  peak  from  the 
back  of  the  TRWE  because  of  the  combined  venting  action  of  the  DSV  and  the  more  open 
P-RWE.  By  the  final  point  in  time  shown  on  this  plot,  186.4  ms,  the  net  overpressure  is  at 
253.1  kPa  (36.7  psi),  still  decreasing,  and  relatively  close  to  the  undisturbed  peak  incident 
shock  overpressure  of  241.3  kPa. 

b.  Force  on  the  Primary  RWE 

Figure  7a  shows  the  axial  forces  on  the  P-RWE  for  the  low  vent  problem,  including  the 
front  face  force,  back  face  force,  and  net  force.  Except  for  the  scaling  factor  resulting  from 
the  multiplication  by  the  area,  the  curves  are  similar  to  those  for  the  average  overpressure. 
The  first  peak  net  force  on  the  P-RWE  is  86.8  MN  (19,500,000  lbf)  at  85.4  ms,  and  the 
second  peak  force  is  92.8  MN  (20,900,000  lbf)  at  108.4  ms.  This  may  be  impractical  to 
design  to,  especially  if  a  safety  factor  of  at  least  2.0  is  applied. 

Figure  7b  shows  a  corresponding  set  of  curves  for  axial  forces  on  the  P-RWE  for  the 
normal  vent. problem.  Here,  the  peak  forces  axe  not  only  substantially  less  than  for  the  low 
vent  problem,  but  the  relief  action  is  also  more  rapid,  as  pointed  out  above.  The  peak  net 
force  on  the  P-RWE  for  the  normal  vent  problem  is  24.7  MN  (5,550,000  lbf),  occurring  at 
85.4  ms.  This  peak  net  force  is  26.6  percent  of  that  for  the  low  vent  case,  due  to  the  combined 
effects  of  having  a  lower  peak  overpressure  and  a  smaller  presented  area,  and  should  be  a 
reasonably  solvable  design  problem. 

2.  Interior  RWE 
a.  Overpressure  on  the  Interior  RWE 

The  average  overpressure  on  the  I-RWE  for  the  low  vent  problem  is  shown  in  Figure  8a, 
illustrating  the  front  face  and  back  face  averages,  and  the  net  value  obtained  by  subtracting 
the  back  face  from  the  front  face.  The  net  overpressure  curve  shows  several  interesting 
features.  There  is  a  peak  value  of  658.8  kPa  (95.6  psi)  at  75.9  ms,  which  results  from  the 
partial  reflection  of  the  incident  shock  from  the  front  face  of  the  I-RWE,  a  rapid  decrease  to 
344.7  kPa  (50.0  psi)  at  78.4  ms,  a  slow  decrease  to  295.8  kPa  (42.9  psi)  at  94.4  ms,  a  rapid 
decrease  to  a  net  negative  value  [in  a  directional  sense]  of  [-]225.1  kPa  ([-]32.6  psi)  at  98.4  ms, 
and  a  small  positive  value  of  7.4  kPa  (1.1  psi)  by  186.4  ms.  The  individual  curves  for  the 
front  and  back  face  show  that  this  behavior  is  caused  by  a  strong  pressure  wave  moving 
upstream.  (Note  the  sudden  rise  in  back  face  average  overpressure  starting  at  94.4  ms.) 
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The  average  overpressure  on  the  I-RWE  for  the  normal  vent  problem  is  shown  in  Fig¬ 
ure  8b.  It  is  immediately  apparent  that  the  curves  for  the  normal  vent  problem  are  distinctly 
different  from  those  for  the  low  vent  problem.  The  initial  peak  overpressure  for  normal  vent¬ 
ing  is  585.9  kPa  (85.0  psi),  which  also  occurs  at  75.9  ms,  but  is  less  than  that  for  low  venting 
because  for  normal  venting  the  incident  shock  must  pass  by  the  open  USV  prior  to  reaching 
the  front  face  of  the  I-RWE.  The  incident  shock  is  thereby  weakened  by  expansion  waves 
moving  .radially  inward  from  the  USV;  the  overpressure  from  the  partial  reflection  of  the 
incident  shock  on  the  I-RWE  front  face  is  also  relieved  more  quickly  for  the  same  reason.  A 
near-plateau  around  300  kPa  (43.5  psi)  is  reached  in  the  net  overpressure  for  the  normal  vent 
problem,  lasting  until  about  98.4  ms.  This  is  followed  by  a  rather  rapid  decrease  to  129.0  kPa 
(18.7  psi),  then  a  slight  increase  and  subsequent  slow  decrease  to  118.3  kPa  (17.2  psi)  by  the 
last  point  in  time  of  186.4  ms.  This  is  caused  by  a  sequence  of  reflected  compression  waves 
and  expansion  waves  from  the  P-RWE,  and  expansion  waves  from  the  USV  and  DSV. 


b.  Force  on  the  Interior  RWE 

The  various  curves  for  axial  forces  on  the  I-RWE  for  the  low  vent  problem  are  shown 
in  Figure  9a,  including  front  face  force,  back  face  force,  and  net  force.  As  noted  previously 
in  the  discussion  for  the  P-RWE,  these  curves  are  simply  the  average  overpressure  curves 
(see  Figure  8a)  scaled  by  area.  The  most  interesting  things  that  this  set  of  curves  highlights 
are  the  reversals  of  force  on  the  I-RWE  as.  time  progresses.  The  peak  net  positive  force 
(acting  toward  the  P-RWE)  is  21.9  MN  (4,920,000  lbf)  at  75.9  ms,  followed  by  a  later  force 
■reversal  to  a  negative  (acting  toward  the  drivers)  net  force  of  7.48  MN  (1,680,000  lbf)  at 
98.4‘  ms  caused  by  the  strong  wave  reflection  from  the  P-RWE.  While  these  forces  do  not  in 
themselves  pose  a  paticularly  difficult  design  problem,  care  would  be  needed  to  design  for  a 
force  reversal. 

Figure  9b  shows  a  corresponding  set  of  curves  for  axial  forces  on  the  I-RWE  for  the 
normal  vent  problem.  There  is  no  force  reversal  here  because  there  is  no  strong  reflected 
wave  from  the  P-RWE.  The  peak  reflected  net  force  on  the  I-RWE  is  19.5  MN  (4,380,000  lbf) 
at  75.9  ms,  caused  by  the  initial  interaction  with  the  incident  shock  wave.  This  peak  then 
decays  in  the  manner  described  for  the  overpressure  curve  for  this  case.  These  forces  do  not 
appear  to  pose  any  particularly  difficult  design  problems,  except  for  the  need  to  properly 
design  the  anchors  for  the  I-RWE  in  the  expansion  tunnel  wall  to  handle  both  the  forces  and 
the  moments  that  they  generate. 

IV.  Mass  Flux  in  the  LB/TS 

1.  Net  Mass  Flux  in  the  LB/TS  as  a  System 

A  discussion  of  the  mass  flow  in  the  LB/TS  is  warranted  here  because  it  helps  to 
illustrate  the  extent  to  which  the  flows  in  the  normal  vent  and  low  vent  computations  are 
different  from  one  another.  For  the  immediate  purpose  of  this  discussion,  it  is  useful  to 
consider  the  LB/TS,  axially  from  the  test  station  (at  20.0  m)  to  the  end  of  the  P-RWE,  and 
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radially  from  the  axis  of  symmetry  to  the  side  wall,  as  a  control  volume.  Figure  10a  shows 
a  set  of  mass  flux  curves  for  the  low  vent  problem.  In  the  discussion  that  follows,  positive 
mass  flux  is  flow  into  the  control  volume  and  negative  mass  flux  is  flow  out  of  the  control 
volume.  The  mass  flux  into  the  control  volume  in  this  case  is  the  mass  flow  across  the  test 
station,  which  is  always  positive  and  decreases  monotonically  with  time.  The  negative  phase 
that  would  occur  in  a  real  blast  wave  is  ignored  in  this  study.  The  low  vent  problem  has 
no  side  venting,  so  the  mass  flow  out  is  only  through  the  restricted  20  percent  opening  in 
the  P-RWE:  The  start  of  mass  flow  through  the  P-RWE  can  be  seen  in  the  curve  labeled 
“OUT”  in  Figure  10a  by  a  change  to  negative  mass  flow  from  no  flow  at  83.4  ms.  The  net 
mass  .flow  curve  shows  immediately  why  the  overpressure  and  force  curves  behaved  as  they 
did.  During  all  of  the  time  simulated  so  fax,  there  is  a  large  mass  accumulation  due  to  the 
positive  net  mass  flow  into  this  control  volume,  specifically  in  the  RWE  end  of  the  LB/TS. 
The  peak  mass  flow  rate  into  the  control  volume  occurs  upon  the  passage  of  the  incident 
shock  by  the  test  station.  The  peak  mass  flow  rate  shown  here  in  Figure  10a,  124.6  .Mg/s 
(275,000  lbm/s),  is  less  than  the  mass  flow  rate  immediately  behind  the  shock  as  it  passes  the 
test  station  because  both  computations  were  started  with  the  shock  well  past  the  test  station, 
just  upstream  from  the  closed  USV.  A  high-pressure,  high-derisity  region  is'  still  expanding 
in  the  upstream  direction  at  186.4  ms,  the  current  ending  time  of  the  computation,  with  the 
mass  still  accumulating  at  the  rate  of  82.5  Mg/s  (182,000  lbm/s). 

Figure  10b  shows  a  similar  set  of  curves  for  mass  flux  for  the  normal  vent  case,  showing 
mass  flux  in,  mass  flux  out,  and  net  mass  flux.  A  comparison  of  these  curves  with  those  in 
Figure  10a  shows  why  the  flows  in  the  two  problems  are  so  different  from  one  another.  The 
normal  vent  problem  has  only,  a  relatively  brief  mass  accumulation  period  because  the  blast 
wave  almost  immediately  passes  by  the  open  USV,  and  then  quickly  reaches  the  open  DSV, 
and  finally  the  more  open  P-RWE.  (The  individual  mass  flows  through  the  test  station,  side 
vents,  and  P-RWE  are  discussed  in  Subsection  IV. 3  below.)  The  net  mass > flux  becomes 
negative  at  85.4  ms,  very  soon  after  the  incident  shock  arrives  a,t  the  P-RWE.  It  stays 
negative  for  the  remainder  of  the  computation  time. 

2.  Mass  Flux  through  the  Interior  RWE 

Figure  11a  shows  the  mass  flux  through  the  I- RWE  for  the  low  vent  problem,  and  gives 
another  good  illustration  of  the  action  of  the  strong  reflected  wave  which  is  moving  upstream. 
The  peak  mass  flow  through  the  I-RWE  is  121.8  Mg/s  (269,000  lbm/s),  occurring  at- 75.9  ms, 
which  is  coincident  with  the  arrival  of  the  incident  shock  at  the  I-RWE.  This  is  not  equal  to 
the  peak  mass  flux  of  the  blast  wave  itself  because  of  the  20  percent  blockage  of  the  I-RWE. 
The  slow  decrease  of  the  mass  flux  with  time  is  relatively  constant  until  98.4  ms,,  and  is  due 
primarily  to  a  combination  of  the  decay  rate  of  the  blast  wave  and  the  partial  reflection  of 
the  incident  shock  from  the  I-RWE.  After  this  time,  the  strong  reflected  wave  (see  Figure  8a) 
from  the  P-RWE  arrives  at  the  back  face  of  the  I-RWE,  decreasing  the  mass  flow  rate  from 
115.0  Mg/s  (254,000  lbm/s)  at  98.4  ms  to  44.0  Mg/s  . (97, 000  lbm/s)  at  104.4  ms,  and  to 
23.6  Mg/s  (52,000  lbm/s)  by  the  ending  time  of  186.4  ms. 

The  curve  for  mass  flow  through  the  I-RWE  for  the  normal  vent  problem,  shown  in 
Figure  lib,  is  markedly  different  from  that  for  the  low  vent  problem.  It  has  a  peak  value 
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of  113.2  Mg/s  (250,000  lbm/s)  at  76.4  ms,  somewhat  below  that  for  the  low  vent  problem 
because  of  the.  action  of  the  USV.  The  curve  then  shows  a  relatively  smooth  monotonic 
decline  with  time,  with  only  a  slight  perturbation  around  110.  ms. 

3.  Individual  Mass  Fluxes  for  the  Normal  Vent  Computation 

Figure  12  shows  the  mass  fluxes  for  each  individual  flow  area  for  the  normal  vent 
problem.  The  curve  labeled  “USV”  represents  the  mass  flux  through  the  USV;  its  greatest 
mass  flux  out  is  35.3  Mg/s  (77,800  lbm/s).  The  curve  labeled  “DSV”  represents  the  mass 
flux  through  the  DSV;  its  greatest  mass  flux  out  is  14.6  Mg/s  (32,200  lbm/s).  The  peak 
mass  flux  out  for  the  P-RWE  is  89.3  Mg/s  (197,000  lbm/s)  at  86.4  ms  corresponding  to  the 
incident  shock  arrival  at  the  P-RWE.  The  mass  flux  at  the  test  station  is  also  shown  for 
reference. 


V.  Overpressure  on  LB/TS  Side  Wall  Sections 

1.  Side  Wall  between  the  Downstream  Side  Vent  and  the  Primary  RWE 

Figure  13a  shows  the  average  overpressure  on  the  LB/TS  side  wall  between  the  closed 
DSV  and  the  P-RWE  for  the  low  vent  problem.  The  solid  curve  shows  the  net  pressure  on 
the  wall.  The  first  peak  of  658.3  kPa  (95.5  psi)  at  91.4  mis  is  followed  by  a  second  peak  of 
677.3  kPa  (98.2  psi)  at  108.4  ms,  and  then  a  gradual  decrease  from  that  point.  The  average 
pressure  is  still  quite  high  at  513.3  kPa  (74.4  psi)  at  the  ending  time  of  186.4  ms,  still  over 
two  times  the  incident  shock  overpressure.  This  is  due  to  the  combined  effects  of  the  pressure 
rise  from  the  reflection  of  the  incident  shock  wave  from  the  P-RWE  and  the  limited  venting. 

Figure  13b  shows  the  average  overpressure  for  the  same  wall  section  for  the  normal  vent 
problem.  Here,  there  is  a  single  peak  overpressure  of  347.6  kPa  (50.4  psi)  at  92.4  ms,  48.7 
percent  less  than  the  highest  peak  for  the  low  vent  problem.  The  decay  from  the  peak  for 
the  normal  vent  problem  is  more  rapid  than  that  for  the  low  vent  problem.  The  last  value 
for  overpressure  at  186.4  ms  is  198.0  kPa  (28.7  psi). 


2.  Side  Wall  between  the  Interior  RWE  and  the  Downstream  Side.  Vent 

Figure  l4a  shows  the  average  overpressure  on  the  LB/TS  side  wall  between  the  I- RWE 
and  the  upstream  end  of  the  closed  DSV  for  the  low  vent  problem;  The  first  peak  in  the  net 
pressure  curve,  214.6  kPa  (31.3  psi)  at  80.4  ms,  is  due  to  the  passage  of  the  incident  shock 
front  after  it  passed  through  the  I- RWE.  After  a  brief  decline  in  pressure,  there  is  a  rapid 
rise  due  to  the.  movement  upstream  of  the  reflected  wave  from  the  P-RWE,  reaching  a  peak 
of  649.1  kPa  (94.1  psi)  at  116.4  ms,  and  then  decreasing  from  that  point. 

Figure  14b  shows  a  similar  set  of  overpressure  curves  for  the  normal  vent  problem.  A 
comparison  between  Figures  14a  and  14b  shows  the  great  difference  in  the  net  side  wall 
pressure  for  the  normal  vent  problem  as  opposed  to  the  low  vent  problem.  The  first  peak 
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net  pressure  is  177.4  kPa  (25.7  psi)  at  80.4  ms,  followed  by  a  subsequent  peak  of  241.2  kPa 
(35.0  psi)  at  106.4  ms,  and  then  a  gradual  decrease  to  171.4  kPa  (24.9  psi)  by  186.4  ms. 

3;  Side  Wall  between  the  Upstream  Side  Vent  and  the  Interior  RWE 

Figure  15a  shows  a  set  of  curves  for  average  overpressure  on  the  LB/TS  side  wall 
between  the  closed  USV  and  the  I-RWE  for  the  low  vent  problem.  After  an  initial  peak  of 
376.7  kPa  (54.6  psi)  at  82.4  ms  due  to  the  passage  of  the  incident  shock  and  its  subsequent 
partial  reflection  from  the  front  face  of  the  I-RWE,  the  net  pressure  reaches  a  maximum  of 
660:0  kPa  (95.7  psi)  at  119.4  ms.  This  sustained  second  high  pressure  is  caused  by  the  tube 
filling  process  discussed  earlier.  The  net  pressure  is  512.3  kPa  (74.3  psi)  at  the  186.4  ms 
ending  time. 

Figure  15b  shows  a  corresponding  set  of  overpressure  curves  for  the  normal  vent  prob¬ 
lem.  The  initial  peak  net  pressure  is  368.4  kPa  (53.4  psi)  at  83.4  ms,  relatively  close  to  the 
initial  peak  for  the  low  vent  problem.  There  is  an  initial  decrease,  a  second  peak  of  354.2  kPa 
(51.4  psi)  at  114.4  ms,  and  then  a  decrease  to  265.3  kPa  (38.5  psi)  by  186.4  ms,  which  is 
considerably  lower  than  the  comparable  net  pressure  for  the  low  vent  problem.  The  initial 
decreases  are  due  primarily  to  the  relieving  action  of  the  USV,  and  the  subsequent  decreases 
are  due  to  the  combined  effects  of  the  larger  vent  areas  relative  to  the  low  vent  case. 

4.  Side  Wall  from  the  Start  of  the  Thick  Wall.  Section  to  the  Upstream  Side 

Vent 

Figure  16a  shows  a  set  of  curves  for  the  overpressure  on  the  LB/TS  side  wall  from 
the  start  of  the  -thick  wall  section  to  the  start  of  the  closed  USV  for  the  low  vent  problem. 
The  initial  decline  in  the  net  pressure  from  225.1  kPa  (32.6  psi)  at  70.6  ms  to  223.4  kPa 
(32.4  psi)  at  89.4  ms  is  due  to  the  normal  decay  of  the  blast  wave.  There  is  a  subsequent  set 
of  two  rises,  leading  to  a  peak  of  590.9  kPa  (85.7  psi)  at  146.4  ms.  The  first  rise  is  from  the 
reflection  of  the  incident  shock  from  the  I-RWE,  and  the  second  rise  is  from  its  reflection 
from  the  P-RWE.  There  is  then  a  decline  to  515.2  kPa  (74.7  psi)  by  186.4  ms  due  to  the 
venting  through  the  P-RWE. 

Figure  16b  shows  a  corresponding  set  of  curves  for  the  normal  vent  problem.  The 
net  pressure  curve  shows  only  minor  disturbances  starting  at  approximately  110  ms;  it  is 
otherwise  similar  to  what  the  pressures  would  be  in  a  very  long  LB/TS  with  no  internal  area 
changes  or  reflecting  surfaces.  The  initial  peak  average  net  pressure  is  233.6  kPa  (33.9  psi) 
at  71.7  ms,  followed  by  a  monotonic  decrease  to  188.3  kPa  (27.3  psi)  by  186.4  ms.  Thus,  it 
appears  that  the.  combined  settings  of  open,  areas  for  the  I-RWE,  P-RWE,  USV,  and  DSV 
allow  the  passage  of  the  blast  wave  through  the  LB/TS  with  a  minimum  of  disturbance. 
This  represents  one  possible  approach,  but  not  necessarily  the  only  means  of  allowing  the 
blast  wave  to  exit  the  expansion  section  without  generating  perturbations  in  the  wave  that 
could  propagate  upstream  to  the  test  station. 
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5.  Entire  Thick  Side  Wall  Section,  Excluding  Open  or  Closed  Side  Vents 

Figure  17a  shows  the  average  overpressure  curves  for  the  entire  length  of  the  thick 
LB/TS  side  wall  section,  excluding  the  closed  side  vents,  for  the  low  vent  problem.  The 
curve  for  net  average  pressure  shows  a  long  steady  rise  from  132.3  kPa  (19.2  psi)  at  70.6  ms 
to  a  peak  of  593.3  kPa  (86.1  psi)  at  141.4  ms,  then  a  decrease  to  513.9  kPa  (74.5  psi)  by 
186.4  ms.  Thus,  there  is  a  relatively  large  average  overpressure,  well  above  the  peak  incident 
shock  overpressure  of  241.3  kPa  (35.0  psi),  in  this  region  for  an  extended  period. 

Figure  17b  shows  a  similar  set  of  curves  for  the  same  side  wall  section,  excluding  open 
side  vents,  for  the  normal  vent  problem.  Here,  the  net  pressure  curve  shows  the  effects  of 
the  open  USV  and  DSV,  and  the  more  open  P-RWE,  by  showing  a  great  relative  reduction 
in  the  average  overpressure  as  compared  with  that  for  the  low  vent  problem.  The  peak  net 
pressure  for  the  normal  vent  problem  is  218.7  kPa  (31.7  psi)  at  91.4  ms,  with  a  value  of 
185.8  kPa  (26.9  psi)  at  the  186.4  ms  ending  time. 


VI.  Test  Station 


Figure  18a  shows  the  average  overpressure  in  the  normal  cross-section  of  the  LB/TS  at 
the  test  station  location  for  the  low  vent  problem  (at  the  20.0  m  point  in  this  computation). 
The  first  value  of  189.3  kPa  (27.5  psi)  at  70.6  ms  is  below  the  peak  incident  shock  overpressure 
of  241.3  kPa  (35.0  psi)  of  the  decaying  blast  wave  because,  at  the  problem  initiation,  this 
point  is  41.9575  m  behind  the  shock  front,  so  it  is  already  exhibiting  a  considerable  amount 
of  decay  in  its  flow  parameters.  Even  though  there  is  a  strong  reflected  shock  still  traveling 
toward  the  test  station  at  the  ending  time  of  186.4  ms,  there  is  as  yet  no  disturbance  from 
it  evident  at  the  test  station.  The  overpressure  at  the  test  station  at  this  time  is  137.1  kPa 
(19.9  psi).  . 

Figure  18b  shows  the  average  overpressure  in  the  same  test  station  cross-section  for 
the  normal  vent  problem.  The  curve  is  identical  to  that  for  the  low  vent  problem  shown  in 
Figure  18a.  Two  reasonable  conclusions  that  can  be  tentatively  drawn  from  this  are  that 
neither  a  compression  wave  (in  the  low  vent  problem)  nor  an  expansion  wave  (in  the  normal 
vent  problem)  from  the  exit  region  of  the  LB/TS  has  yet  reached  the  test  station. 


VII.  Conclusion 


The  results  and  observations  discussed  here  are  to  be  considered  preliminary.  It  does 
appear  that  the  combined  RWE  and  side  vent  design  analyzed  here  is  efficient  in  its  reduction 
of  reflected  compression  and  shock  waves  from  the  RWE’s  while  minimizing  the  generation  of 
excessively  strong  expansion  waves.  The  loadings  generated  on  the  RWE’s  and  the  LB/TS 
side  wall  appear  to  be  within  ranges  that  allow  practical  designs,  assuming  that. all  vent 
areas  are  at  theit  proper  settings. 

The  low  vent  problem,  embodying  a  worst-case  assumption  of  the  possible  settings  of 
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the  areas  of  the  USV,  DSV,  and  P-RWE,  combined  with  the  strongest  design  blast  wave, 
appears  to  generate  blast  loads  on  the  RWE’s  that  are  beyond  the  reach  of  a  reasonable, 
practical  design.  The  loads  on  the  LB/TS  side  wall  may  also  be  prohibitively  large.  It  seems 
reasonable  to  design  the  LB/TS  with  a  set  of  interlocks  with  logic  that  requires  failure  to  a 
safe  state  such  that  the  drivers  cannot  be  pressurized,  and  the  diaphragms  (or  throat  valves) 
cannot  be  burst  (opened)  unless  the  RWE(’s)  and  side  vents  (if  any)  are  in  their  proper 
positions.  It  would  also  be  prudent  to  have  a  design  for  the  RWE(’s)  and  any  side  vents 
such  that  they  would  fail  to  a  fully  open  position  to  further  assure  the  saving  of  the  LB/TS 
in  the  event  of  an  off-design  firing. 
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ELIMINATOR 

Figure  la.  A  Simplified  Schematic  of  the  Proposed  US  LB/TS.  (Courtesy  of  R.  Pearson.) 


Figure  lb.  Proposed  LB/TS  Cross-Section.  (Courtesy  of  R.  Pearson.) 
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Figure  2.  Axial  Section  through  the  Expansion  Tunnel  (Section  A-A  in  Figure  3).  (Courtesy 
of  R.  Pearson.). 
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Figure  3.  Exit  Plane  of  the  Primary  Rarefaction  Wave  Eliminator  (Section  D-D  in  Figure  2) 
(Courtesy  of  R.  Pearson.)  &  1 
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Figure  4.  Two  Candidate  Shapes  of  the  Vanes  for  the  Primary  Rarefaction  Wave  Eliminator. 
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Figure  5.  Tunnel  Cross  Section  Looking  Upstream  at  the  Interior  Rarefaction  Wave  Elim 
inator  (Section  C-C  in  Figure  2).  (Courtesy  of  R.  Pearson.) 
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Figure  7b.  Axial.  Forces  on  the  P-RWE,  Normal  Vent. 
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Figure  8b.  Average  Overpressure  on  the  I-R.WE,  Normal  Vent. 
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Figure  9b.  Axial  Forces  on  the  I-RWE,  Normal  Vent. 
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Figure  11a.  Mass  Flux  through  the  I-RWE,  Low  Vent. 
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Figure  lib.  Mass  Flux  through  the  I-RWE,  Normal  Vent. 
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Figure  13a.  Average  Overpressure  on  the  LB/TS  Side  Wall  between  the  DSV  and  the 
P-RWE,  Low  Vent. 
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Figure  13b.  Average  Overpressure  on  the  LB/TS  Side  Wall  between  the  DSV  and  the 
P-RWE,  Normal  Vent. 
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Figure  14a.  Average  Overpressure  on  the  LB/TS  Side  Wall  between  the  I-RWE  and  the 
DSV,  Low  Vent. 
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Figure  14b.  Average  Overpressure  on  the  LB/TS  Side  Wall  between  the  I-RWE  and  the 
DSV,  Normal  Vent. 
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Figure  15a.  Average  Overpressure  on  the  LB/TS  Side  Wall  between  the  USV  and  the 
I-RWE,  Low  Vent. 
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Figure  15b.  Average  Overpressure' on  the  LB/TS  Side  Wall  between  the  USV  and  the 
I-RWE,  Normal  Vent. 
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Figure  16a.  Average  Overpressure  on  the  LB/TS  Side  Wall  from  the  Start  of  the  Thick 
Wall  Section  to  the  USV,  Low  Vent. 
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Figure  16b.  Average  Overpressure  on  the  LB/TS  .Side  Wall  from  the  Start  bf  the  Thick 
Wall  Section  to  the  USV,  Normal  Vent. 
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.Figure  17a.  Average  Overpressure  on  the  LB/TS  Thick  Side  Wall  Section,  Excluding 
Closed  Side  Vent  Areas,  Low  Vent. 
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Figure  18a.  Average  Overpressure  across  the  LB/TS  Cross-Sectional  Area  at  the  Test 
Station,  Low  Vent. 
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Figure  18b.  Average  Overpressure  across 'the  LB/TS  Cross-Sectional  Area  at  the  Test 
Station,  Normal  Vent. 
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